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ABSTRACT 


Travel  times  of  PcP  and  P  phases  from  nuclear  and  high  explosive 
sources  are  interpreted  in  terms  of  variations  in  the  radius  of  the 
earth's  outer  core.  Travel-time  residuals  from  the  Taggart-Engdahl 
[1968]  tables  are  determined  for  the  PcP-P  time  interval  and  for  PcP 
travel  times,  corrected  for  station  and  source  terms.  The  travel  times 
are  also  corrected  for  elevation  of  source  and  station  and  ellipticity, 
excluding  the  core  ellipticity  term.  This  interpretation  favors  a  core 
which  is  slightly  larger  than  that  of  the  reference  model  of  Taggart  and 
Engdahl  [1968]  and  which  has  less  ellipticity  than  estimated  by  Bullen 
[1936],  Other  interpretations  in  terms  of  lateral  variation  in  mantle 
velocity  are  also  possible. 

Because  all  of  the  events  and  most  of  the  stations  whose  data  are 
used  lie  in  the  northern  hemisphere,  the  description  of  core  shape  is  in¬ 
complete  and  is  representative  only  of  regions  of  the  core  surface  in  the 
vicinity  of  the  PcP  reflection  points. 

The  distribution  of  data  does  not  permit  a  conventional  spherical 
harmonic  analysis  of  the  inferred  variation  in  core  radius.  A  modified 
spherical  harmonic  analysis  method  is  devised  to  smooth  the  data  and  to 
estimate  the  shape  of  the  core.  The  variation  of  core  radius  determined 
from  this  representation  is  approximately  ±10  kilometers  on  the  basis  of 
terms  to  degree  and  order  5.  The  standard  deviation  of  the  data  is 
approximately  ±20  kilometers,  indicating  the  presence  of  variations  of 
higher  degree  and  order,  many  of  which  are  probably  due  to  factors 
other  than  undulation  of  the  mantle-core  boundary.  Variations  in  geoid 
height  resulting  from  the  inferred  variation  in  core  radius  do  not 
correlate  with,  but  are  of  the  same  order  of  magnitude  as,  those  deter¬ 
mined  from  satellite  observations. 

Analysis  of  a  limited  number  of  short-period  PcP  amplitudes  shows 
no  significant  variation  of  the  PcP  amplitude/period  ratio  with  epicentral 


vii 


distance  and  an  estimated  scatter  of  data  on  the  same  order  as  that  found 
for  P  in  the  distance  range  70°  to  90°.  The  PcP  phase  is  generally 
shorter  in  apparent  period  than  the  P  phase.  This  relationship  is  re¬ 
versed,  however,  for  arrivals  from  the  Novaya  Zemlya  event  of  October 
27,  1966. 

The  PcP  amplitude  and  travel- time  variance  can  be  explained,  for 
the  most  part,  in  terms  of  lateral  variations  in  thickness  and  distri¬ 
bution  of  crustal  and  upper  mantle  velocity  zones  consistent  with  the 
model  proposed  by  Nuttli  and  Bolt  [1969]  for  northern  California. 
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Introduction 


The  purpose  of  this  study  is  to  examine  the  variation  in  PcP  travel 
times  and  the  PcP-P  time  interval  from  zero-focus  sources  and  to  estimate 
on  this  basis  the  shape  of  the  core.  Because  of  the  distribution  of 
zero-focus  sources,  the  results  apply  only  to  the  northern  hemisphere. 

The  contributions  of  crust  and  upper  mantle  velocity  variations  to  the 
PcP  and  PcP-P  scatter  are  assessed,  and  appropriate  corrections  applied 
whenever  possible. 

Background 

A.  Vogel  [1961  |  considered  the  travel  time  dispersions  of  the  PcP, 
ScS,  PcS,  and  ScP  waves  as  consequences  of  a  non-uniform  reflection  depth. 
Using  data  from  earthquake  sources,  he  interpreted  travel  time  anomalies 
as  corresponding  to  variations  of  more  than  +200  kilometers  in  the  depth 
to  the  core-mantle  boundary.  The  geographical  pattern  of  core  depths  was 
found  to  differ  for  the  periods  1948-1954,  and  1930-1936.  He  also  found 
correlation  of  the  core  depth  anomalies  with  those  of  the  earth's  gravity 
and  magnetic  fields.  Egyed  [1964]  accounted  for  the  amplitude  of  the  un¬ 
dulations  of  Izsak's  [1964]  satellite  geoid  in  terms  of  the  density  con¬ 
trast  of  a  sphere  of  200  km  diameter  situated  at  the  core-mantle  boundary. 

The  present  study  differs  from  that  of  Vogel  in  several  important 

ways : 

1.  Only  surface-focus  (nuclear  or  high  explosive)  events 

are  used. 

2.  In  each  case  the  zero  time  is  either  known  or  is  assumed 
by  adjusting  the  computed  zero  times  based  on  the  Jeffreys-Bullen  sur¬ 
face-focus  times  by  a  regional  correction.  The  adjusted  times  indicate 
that  the  shots  were  fired  on  or  within  a  fraction  of  a  second  of  an 
exact  minute;  hence  in  most  cases  the  unknown  origin  times  are  assumed 
to  be  exactly  on  the  minute. 

3.  Epicenter  locations  used  in  this  study  are  either  known 
a  priori  or  have  been  determined  with  an  accuracy  unattainable  a  few 
years  ago. 

4.  Depth  phases  are  not  present,  eliminating  a  number  of 
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arrivals  which,  at  certain  distances,  could  be  confused  with  core  phases. 

5.  Only  PcP  and  P  travel  times  from  short-period  instruments 

are  used. 

PcP  Observations  from  Zero-Focus  Sources 

The  data  used  in  this  study  represent  all  available  PcP  readings 
from  zero-focus  events  that  are  known  to  the  author.  To  avoid  the  redun¬ 
dancy  inherent  in  using  all  available  data  from  NTS  sources,  only  the 
Bilby  and  Faultless  data  are  included.  The  source  regions  and  events 
used  are  described  in  Table  1  and  shown  in  Figure  1.  Data  are  included 
from  many  recording  sites,  including  U.S.C.&.G.S. *  U.S.G.S.,  AFTAC 
stations,  LASA,  and  a  number  of  independent  university  or  foreign  govern¬ 
ment-operated  stations.  Epicentral  distances  for  the  PcP  observations 
range  from  15°  to  85°.  Published,  bulletin,  or  listed  arrivals  identified 
as  PcP  have  been  accepted  as  such  and  included  in  the  analysis  unless 
inspection  of  the  record  by  the  author  or  others  has  shown  this  to  be  in¬ 
correct.  In  some  cases  unidentified  arrivals  or  misidentified  arrivals 
are  included  when  the  travel  time  indicates  the  phase  should  be  PcP. 

Data  from  western  Pacific  sources  were  obtained  from  Kogan  [1960], 
Carder  [1964],  Carder  [1968^,  and  Doyle  and  Webb  [1958],  Buchbinder 
[1965]  and  Jordan  [1968]  supplied  the  Bilby  and  Faultless  data  respective¬ 
ly.  Longshot  data  were  from  Carder,  et  a!.  [1967],  Buchbinder  [1968a], 
Taggart  and  Engdahl  [1968],  Clark  [1966a],  Bufe  and  Willis  [1969],  and 
others.  CHASE  III  and  IV  arrivals  are  from  Buchbinder  [1969],  Clark 
[1965],  Clark  [1966b],  and  from  the  ESSA  Earthquake  Data  Reports.  The 
travel  times  for  the  Sahara  event  were  read  by  the  author  from  records  of 
the  ESSA  worldwide  net.  Semipalatinsk  and  Novaya  Zemlya  data  are  from 
Jordan  [1968],  from  B.C.I.S.,  ESSA,  and  LASA  bulletins,  and  from  the 
author's  readings.  All  PcP  travel  times  used  are  tabulated  in  Appendix 
1,  where  the  sources  of  data  are  referenced  in  detail.  Readings  made  or 
verified  by  the  author  are  indicated  by  an  asterisk.  Arrivals  determined 
by  the  author  have  been  subjected  to  several  criteria,  which  follow. 

1.  Arrival  should  be  distinct  from  surrounding  microseismic 
background  and  P  coda. 

2.  Generally  PcP  will  be  equal  to  or  higher  than  P  in  frequency 
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Table  1.  Description  of  Surface-Focus  Events 
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FIGURE  I.  DISTRIBUTION  OF  ZERO-FOCUS  SOURCES. 


at  distances  greater  than  30°.  This  relationship  has  been  noted  by 
Kogan  [1968]  and  is  explained  by  Kanamori  [1967]  in  terms  of  attenuation 
characteristics  in  the  mantle.  Arrivals  from  the  Fovaya  Zemlya  event  of 
10/27/66  are  exceptions  to  this  rule. 

3.  At  epicentral  distances  greater  than  about  30°,  the  wave¬ 
forms  of  PcP  and  P  should  be  similar,  such  that  the  PcP-P  time  interval 
could  be  determined  by  auto-correlation.  Spectral  analyses  by  Kanamori 
fl967]  and  Buchbinder  [1969]  have  shown  this  to  be  true. 

4.  In  many  instances  for  stations  at  the  shorter  epicentral 
distances  both  horizontal  and  vertical  component  seismograms  were  scan¬ 
ned.  The  ratio  of  vertical  to  horizontal  amplitudes  should  be  larger 
for  PcP  than  for  P. 

Since  many  of  the  arrivals  selected  by  the  author  did  not  meet  all 
of  these  criteria  and  the  reliability  of  the  bulletin  data  in  this  re¬ 
spect  is  unknown,  no  attempt  was  made  to  assign  weights  to  readings  on 
this  basis. 

Correction  for  the  Earth's  Ellipticity 

Jeffreys  [1935]  described  the  effect  of  the  earth's  ellipticity  on 
the  travel  times  of  seismic  body  waves.  The  ellipticity  correction  for 
P  or  S  waves  involves  terms  dependent  upon  the  ellipticity  of  the  earth's 
outer  surface  and  the  latitudes  of  source  and  station  and  upon  the  ellip- 
ticities  of  internal  strata  of  equal  velocity  along  the  ray  path.  Bullen 
[1936]  assumed  hydrostatic  conditions  to  apply  within  the  earth  and  cal¬ 
culated  the  variation  of  density  with  depth  on  the  basis  of  the  velocities 
of  P  and  S  waves.  The  ellipticities  of  internal  strata  of  equal  density 
were  then  computed,  ranging  from  0.00337  at  the  earth's  surface  to  0.00256 
at  its  center. 

The  ellipticity  correction  to  the  PcP  travel  time  can  be  considered 

as  four  separate  terms  with  geometry  as  illustrated  in  Figure  2.  Jeffreys 

6  r 

[1935]  has  shown  that,  to  the  first  order  in  —  ,  the  difference  between 
the  actual  travel  time  of  the  P  or  S  ray  between  source  and  station  and 
the  time  between  the  two  corresponding  points  on  the  surface  of  the  stan¬ 
dard  sphere,  is  given  by 
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(T'-T) 


P  or  S 


(1) 


where  p  is  the  ray  parameter  defined  at  any  point  on  the  ray  path  as 


d 0 

p  =  —  —  ,  v  is  the  local  wave  velocity,  and 


(2) 


r  is  the  distance  of  the  point  from  the  center  of  the  earth.  Consider¬ 
ing  the  PcP  ray  path  in  two  segments,  source  to  core  and  core  to  station, 
each  segment  subtending  a  central  angle  A/2,  the  above  expression  becomes 


The  first  two  terms  of  (3)  represent  the  "end  effects"  portion  of 
the  ellipticity  correction.  By  (2),  they  can  be  written 


(4) 


Since  —  is  negative  on  the  down  path  (source  to  core),  (4)  can  be  writ- 
ds 

ten 


r£2S-i6r]  _  2  [ass-i  Sri  „  +  N-i  J 

v  J  station  |_  v  J  mantle-core  [  v  J  source  • 

boundary 


(6) 
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This  is  the  standard  "delay  time"  format  often  used  in  refraction 
studies  and  in  applying  elevation  corrections  to  seismic  travel  times. 

The  second  term  of  (6)  has  been  excluded  from  the  ellipticity  corrections 
applied  in  this  work,  since  it  is  6r  at  the  core-mantle  boundary  which  we 
wish  to  determine,  i.e.,  for  the  purpose  of  computing  ellipticity  correc¬ 
tions,  the  core  is  assumed  to  be  spherical. 

Thus  the  "end  effects"  portion  of  the  ellipticity  correction  is  re¬ 
duced  to  the  source  and  station  terms  of  (6),  which  can  be  simplified  to 

I  cos-  i|  •  (Sr  +  5r  ),  assuming  (7) 

I  v  I,  r  source  station  *  6 

L  Jbase  of 

crust 

similar  conditions  at  opposite  ends  of  the  ray  path. 

The  path  difference  is  assumed  to  occur  in  the  upper  mantle.  The 
factor  cos  i/v  at  the  base  of  the  crust  has  been  evaluated  as  a  function 
of  epicentral  distance  for  PcP  and  is  plotted  in  Figure  3.  Angle  of  in¬ 
cidence  was  determined  from  the  relationship 

sin  1  "  7  ft  (8) 

assuming  a  sub-crustal  velocity  of  8.1  km/sec.  The  values  of  St/60  were 
taken  from  Taggart  and  Engdahl's  [1968]  PcP  travel-time  tables. 

The  Sr  values  in  equation  (7)  are  functions  of  the  latitude  of  the 
source  or  station  and  the  ellipticity  of  the  earth's  surface.  In  general, 
from  the  definition  of  the  standard  sphere 

Sr  =  er  S2  (9) 

where  e  is  the  ellipticity  and  S2  is  the  surface  harmonic,  S2  =  -j  -  cos2<f>, 
where  4>  is  the  geocentric  colatitude.  The  Sr  sum  was  computed  for  each 

source-station  combination  and  used  in  conjunction  with  Figure  3  and  Eq. 
(7)  to  compute  the  "end  effects"  portion  of  the  ellipticity  correction. 

Elevation  corrections  were  applied  in  the  same  manner,  assuming  a  crustal 
velocity  of  6.0  km/sec  and  cos  i/v  versus  epicentral  distance  as  plotted 
in  Figure  4. 
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FIGURE  4.  DELAY  FACTOR,  ELEVATION  CORRECTION. 
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The  integral  terms  of  (3)  represent  the  part  of  the  correction  due 
to  the  ellipticity  of  internal  strata  along  the  ray  path.  Bullen  [1938], 
using  his  [1936]  values  of  ellipticity  of  internal  strata,  has  shown  that 
this  expression  can,  subject  to  errors  not  exceeding  0.1  second,  be  put 
in  the  form  +  y  Sa  +  S^) ,  where  and  are  the  values  of  the 

surface  harmonic  at  the  epicenter  and  observatory  respectively,  and 
S  the  value  at  the  point  of  reflection  on  the  core-mantle  boundary. 

cL  ^ 

Bullen' s  computed  values  of  fj^)  are  P^otte^  as  a  function  of  A  in  Figure 
5. 

Bullen's  [1936]  values  of  core  radius  and  core  ellipticity  are  in¬ 
volved  in  tnis  approximation,  but  because  of  the  integral  nature  of  this 
term,  perturbations  of  the  core  radius  and  ellipticity  of  the  order  en¬ 
countered  in  this  study  would  not  significantly  change  the  correction. 

2 

The  term  (Sn  +  —  S  +  S..)  was  computed  for  each  epicenter  and  station 

U  J  3  -L 

combination  and  used  in  conjunction  with  the  curve  of  Figure  5  to  compute 
the  correction  due  to  the  ellipticity  of  internal  strata.  The  summed 
ellipticity  corrections  were  then  combined  with  elevation  corrections  for 
each  source  and  station  combination,  applied  to  the  observed  travel  times, 
and  the  resulting  corrected  travel  times  compiled  in  Appendix  1.  Care 
was  taken  in  computing  the  ellipticity  correction  to  avoid  the  introduc¬ 
tion  of  systematic  errors  which  would  affect  the  determination  of  the 
core's  shape. 

Determination  of  Travel~Time  Residuals 

The  datum  chosen  for  this  study  is  that  of  Taggart  and  Engdahl  [1968], 
whose  curve  is  based  on  the  model  of  Herrin,  et  al.  [1968]  and  a  core 
depth  determined  from  zero-focus  PcP  travel  times.  Residuals  from  the 
Taggart -Engdahl  tables  were  determined  for  the  corrected  travel  times. 
Another  set  of  residuals  was  determined  by  applying  station  and  source 
corrections.  These  were  applied  in  the  same  manner  used  by  Taggart  and 
Engdahl  [1968],  assuming  nearly  vertical  incidence  and  neglecting  azimu¬ 
thal  terms.  Station  and  source  corrections  with  appropriate  datum  adjust¬ 
ments  were  taken  from  Herrin  and  Taggart  [1968],  Cleary  and  Hales  [1966], 
Carder,  Gordon,  and  Jordan  [1966  ],  estimated  from  Figure  2  of  Toksoz  and 
Arkani  Hamed  [1967],  or  determined  by  the  writer.  It  was  hoped  that  the 
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scatter  due  to  heterogeneity  of  the  crust  and  upper  mantle  would  be  re¬ 
duced  by  the  application  of  these  corrections.  Another  possible  method 
of  accomplishing  this  is  by  the  use  of  the  PcP-P  time  interval,  the 
method  used  by  Vogel  [I960],  The  use  of  PcP-P  hes  many  advantages  and  a 
few  disadvantages.  At  distances  beyond  about  35°  the  paths  of  P  and  PcP 
through  the  crust  and  upper  mantle  are  close  enough  together  to  produce 
near  cancellation  of  station  and  source  effects  for  the  two  phases  in 
most  cases.  Errors  in  station  time  or  origin  time  will  not  affect  PcP-P. 

On  the  negative  side,  variations  in  compressional  wave  velocity  near  the 
"bottom"  of  the  P  ray  path  would  produce  P-wave  travel-time  anomalies  but 
would  have  little  effect  on  PcP,  thus  changing  the  PcP-P  interval.  In 
addition,  while  P  is  normally  a  stronger  phase  than  PcP,  there  are  some 
instances  where  the  onset  of  PcP  is  more  sharply  defined  than  that  of  P. 

Two  such  examples  are  shown  in  Figure  6. 

In  order  to  determine  the  PcP-P  arrival-time  interval,  P-wave  travel 
times  as  read  and  corrected  for  ellipticity  and  elevation,  have  been  tabu¬ 
lated  with  the  PcP  times  and  are  included  in  Appendix  1.  Residuals  from 
the  Taggart-Engdahl  [1968]  PcP-P  table  have  been  calculated,  with  correc¬ 
tions  applied  for  ellipticity  and  elevation. 

The  three  sets  of  residuals  described  above  are  tabulated  in  Appen¬ 
dix  3  and  are  shown  as  functions  of  epicentral  distance  in  Figures  7 
through  12.  Figures  7  and  8  show  PcP  residuals,  corrected  for  ellipticity 
and  elevation  only,  for  U.  S.  events  and  for  foreign  events,  respectively. 
Figures  9  and  10  are  the  same  presentations  with  additional  source  and 
station  corrections  applied  to  the  residuals.  PcP-P  residuals  are  shown 
in  Figure  1.1  for  U.S.  explosive  sources,  and  in  Figure  12  for  foreign  events. 

Estimate  of  Travel-Time  Error 

There  are  several  sources  of  errors  which  will  affect  the  PcP  travel¬ 
time  residual  determinations.  These  errors  are  of  two  types;  those  which 
can  be  estimated,  and  those  which  cannot.  The  first  class  are  listed 
below  in  Tables  2  for  PcP  residuals  and  3  for  PcP-P,  with  estimates  of 
their  magnitudes.  Error  estimates  are  considered  for  two  ranges  of  epi¬ 
central  distance  (<45°  and  >45°)  and  source  areas  are  classified  accord¬ 
ing  to  the  believed  accuracy  of  zero  time  and  epicenter  location. 
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FIGURE  6. 
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FIGURE  7.  PcP  RESIDUA  1*3  FROM  TAGGART-ENGDAHL  TABLES,  UNITED  STATES  EVENTS. 


i 


(SONOD3S)  nwiaissa  d^d 


16 


FIGURE  8.  PcP  RESIDUA  15  FROM  TAGGART-ENGDAHL  TABLES,  FOREIGN  EVENTS 
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FIGURE  9.  PcP  RESIDUALS  AFTER  APPLICATION  OF  SOURCE  AND  STATION  CORRECTIONS,  UNITED  STATES  EVENTS 
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FIGURE  10.  PcP  RESIDUALs  AFTER  APPLICATION  OF  SOURCE  AND  STATION  CORRECTIONS,  FOREIGN 
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FIGURE  12.  PcP-P  RESIDUAIS,  FOREIGN  EVENTS. 


Table  2.  Error  Estimation,  PcP  Residuals 


Typ< 

e  of  Error 

Distance  Range 

Error 

in  Seconds 

for  Event 

U.S. 

French 

Russian 

1. 

Epicenter  Location 

A 

<  45° 

_ 

±  .1 

±  .2 

A 

>  45° 

— 

±  .2 

±  .4 

2. 

Zero  Time 

— 

±  .2 

±  .3 

3. 

Station  Time 

±0.1 

±0.1 

±0.1 

4. 

Error  in  Scaling  Time 

±0.2 

±0.2 

±0.2 

of  PcP  from  Record 

Estimate  of  Mean 

A 

<  45° 

±0.2 

±0.3 

±0.4 

Error  Sum 

A 

>  45° 

±0.2 

±0.4 

±0.5 

Table  3.  PcP-P 

Estimate  of  Error  in  Travel-Time  Residual 


Type  of  Error 

Distance  Range 

Error  in 

Seconds 

for  Event 

U.S. 

French 

Russian 

1. 

Epicenter  Location 

A 

<  45° 

negligible 

±  .3 

±  .5 

A 

>  45° 

negligible 

±  .1 

±  .2 

2. 

Zero  Time 

All 

no  effect  - 

3. 

Station  Time 

All 

no  effect  - 

4. 

Error  in  Scaling  Time 

±0.2 

±0.2 

±0.2 

of  PcP  from  Record 

5. 

Error  in  Scaling  Time 

±0.1 

±0.1 

±0.1 

of  P  from  Record 

Estimate  of  Mean 

A 

<  45° 

±0.2 

±0.4 

±0.5 

Error  Sum 

A 

>  45° 

±0.2 

±0.2 

±0.3 
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Estimates  of  mean  error  sum  are  determined  using  the  method  of  random 
vector  addition: 

Mean  Error  Sum  =  /e, 2+E02+E.2+E, 2+E  * 

V  1  2  3  4  5 

From  Tables  2  and  3  it  is  seen  that  the  sources  of  error  in  the  PcP 
residuals  are  somewhat  different  from  those  involving  PcP-P  residuals. 

For  instance,  the  error  due  to  an  incorrect  epicenter  location  will  be 
more  pronounced  at  short  distances  for  PcP-P  and  at  long  distances  for 
PcP  residuals 

The  second  class  of  error  cannot  be  measured  or  predicted  with  any 
degree  of  certainty,  and  involves  mis-identification  of  arrivals  or  fail¬ 
ure  to  select  the  "first  break"  of  PcP. 

In  regard  to  the  first  possibility,  there  is  often  a  tendency  to 
"seek  out"  a  PcP  by  selecting  the  strongest  arrival  nearest  the  predicted 
PcP  time.  If  the  author's  criteria  for  PcP  discussed  earlier  are  not 
heeded,  and  sometimes  even  if  they  are  heeded,  mis-identification  is  a 
distinct  possinility.  The  primary  problem  here  is  that  the  assumption 
that  a  readable  PcP  phase  must  be  present  somewhere  on  the  seismogram  is 
incorrect.  It  has  been  shown  (see  Carder,  1964)  that  PcP  is  a  "will  of 
the  wisp"  which  may  appear  or  disappear  with  a  slight  shift  of  epicenter 
or  recording  station.  The  wide  scatter  of  PcP  amplitudes  from  Longshot 
(Carder,  et  al. ,  1967  and  Clark,  1966a)  is  further  evidence  of  this  phenom¬ 
enon.  Amplitudes  of  some  of  the  arrivals  whose  travel  times  are  used  here 
are  discussed  in  a  later  section  of  this  paper. 

At  distances  near  45°,  PP  and  PcP  are  sometimes  confused,  although 
examination  of  the  horizontal  component  records  will  usually  resolve  this 
difficulty.  In  a  narrow  distance  range  near  60°  there  appears  to  be  a 
series  of  arrivals  preceding  PcP  and  having  an  apparent  velocity  between 
those  of  first  P  and  PcP.  One  possible  such  phase  is  the  "B"  phase  of 
Gutenberg  and  Richter  [1935]  whose  travel-time  curve  parallels  that  of 
PcP  and  is  9  to  18  seconds  earlier.  The  unidentified  arrival  shown  at  the 
top  of  Figure  13  could  be  taken  as  PcP  were  it  not  for  the  "on  time"  PcP 
arrival  18  seconds  later.  At  distances  beyond  70  degrees,  crustal  rever¬ 
berations  may  produce  pseudo  PcP's.  Longshot  and  Novaya  Zemlya  events 
show  a  complex  P  coda  at  most  distances.  Arrivals  from  separate  events 
are  also  possible  sources  of  a  coincidental  error. 


22 


Event:  Novaya  Zenlya  10/27/66  recorded  at  Edeonton,  Alberta  (EOH),  A  «  53.413° 
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Estimates  of  mean  error  sum  are  determined  using  the  method  of  random 
vector  addition: 

Mean  Error  Sum  =  /e, 2+E02+E„2+E, 2+Ec 2 
V  1  2  3  4  5 

From  Tables  2  and  3  it  is  ssen  that  the  sources  of  error  in  the  PcP 
residuals  are  somewhat  different  from  those  involving  PcP-P  residuals. 

For  instance,  the  error  due  to  an  incorrect  epicenter  location  will  be 
more  pronounced  at  short  distances  for  PcP-P  and  at  long  distances  for 
PcP  residuals 

The  second  class  of  error  cannot  be  measured  or  predicted  with  any 
degree  of  certainty,  and  involves  mis-identification  of  arrivals  or  fail¬ 
ure  to  select  the  "first  break"  of  PcP. 

In  regard  to  the  first  possibility,  there  is  often  a  tendency  to 
"seek  out"  a  PcP  by  selecting  the  strongest  arrival  nearest  the  predicted 
PcP  time.  If  the  author’s  criteria  for  PcP  discussed  earlier  are  not 
heeded,  and  sometimes  even  if  they  are  heeded,  mis-identification  is  a 
distinct  possibility.  The  primary  problem  here  is  that  the  assumption 
that  a  readable  PcP  phase  must  be  present  somewhere  on  the  seismogram  is 
incorrect.  It  has  been  shown  (see  Carder,  1964)  that  PcP  is  a  "will  of 
the  wisp"  which  may  appear  or  disappear  with  a  slight  shift  of  epicenter 
or  recording  station.  The  wide  scatter  of  PcP  amplitudes  from  Longshot 
(Carder,  et  al.  ,  1967  and  Clark,  1966a)  is  further  evidence  of  this  phenom¬ 
enon.  Amplitudes  of  some  of  the  arrivals  whose  travel  times  are  used  here 
are  discussed  in  a  later  section  of  this  paper. 

At  distances  near  45°,  PP  and  PcP  are  sometimes  confused,  although 
examination  of  the  horizontal  component  records  will  usually  resolve  this 
difficulty.  In  a  narrow  distance  range  near  60°  there  appears  to  be  a 
series  of  arrivals  preceding  PcP  and  having  an  apparent  velocity  between 
those  of  first  P  and  PcP.  One  possible  such  phase  is  the  "B"  phase  of 
Outenberg  and  Richter  [1935]  whose  travel-time  curve  parallels  that  of 
PcP  and  is  9  to  18  seconds  earlier.  The  unidentified  arrival  shown  at  the 
top  of  Figure  13  could  be  taken  as  PcP  were  it  not  for  the  "on  time"  PcP 
arrival  18  seconds  later.  At  distances  beyond  70  degrees,  crustal  rever¬ 
berations  may  produce  pseudo  PcP's.  Longshot  and  Novaya  Zemlya  events 
show  a  complex  P  coda  at  most  distances.  Arrivals  from  separate  events 
are  also  possible  sources  of  a  coincidental  error. 
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Event:  Hovaya  Zemlya  10/27/66  recorded  at  Edaonton,  Alberta  (EOM),  6  «  53.413° 
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FIGURE  13.  SURFACE-FOCUS  EVENTS  RECORDED  AT  NORTH  AMERICAN  STATIONS  ON  VERTICAL  SHORT- PERIOD  INSTRUMENTS. 


In  selecting  the  "first  break"  of  PcP,  the  author  generally 
found  the  relationship  of  Buchbinder  [1968a]  to  hold  true,  with  like 
polarity  for  PcP  and  P  at  distances  greater  than  about  30  degrees  and 
opposite  polarity  at  shorter  distances.  "Visual  cross-correlation"  was 
used  to  determine  the  PcP-P  interval,  and  hence  the  PcP  arrival  time,  at 
distances  larger  than  about  30  degrees.  At  shorter  distances  the  first 
distinct  motion  was  selected,  and  was  generally  a  dilation. 

Possible  Interpretations  of  Residuals 

The  travel-time  residuals  described  in  the  previous  section  can  be 
interpreted,  within  the  limitations  of  their  probable  errors,  as  due  to 
several  factors: 

1.  Anomalous  crust  and  upper  mantle  velocities  near  epicenter 
or  station. 

2.  Anomalous  lower  mantle  velocities. 

3.  Variations  in  depth  to  the  core-mantle  boundary. 

4.  Slope  or  "dip"  of  the  core-mantle  interface. 

Factor  1  is  certainly  significant;  station  and  source  anomalies  do 
exist  as  shown  by  Carder,  Gordon,  and  Jordan  [1966],  and  may  be  highly 
directional  as  demonstrated  by  Cleary  and  Hales  [1966],  Cleary  [1967], 
and  Herrin  and  Taggart  [1968],  These  effects  serve  to  increase  the  scat¬ 
ter  of  the  data  and  may  also  produce  systematic  errors.  It  is  believed 
that  such  errors  will  be  reduced  by:  1)  the  application  of  source  and  sta¬ 
tion  corrections  or  by  (2)  the  use  of  PcP-P  times  in  the  analysis. 

Anomalous  lower  mantle  velocities  would  produce  residuals  very  dif¬ 
ficult  to  eliminate  or  distinguish  from  those  due  to  variations  in  core 
depth.  The  velocity  structure  in  the  lower  four  hundred  kilometers  is 
not  known  (Taggart  and  Engdahl,  1968)  with  any  degree  of  certainty  and  may 
vary  with  latitude  and  longitude. 

The  residuals  for  epicentral  distances  less  than  70  degrees  will  be 
interpreted  in  terms  of  variations  in  depth  to  the  core-mantle  boundary. 
Such  data  should  permit  a  determination  of  the  gross  features  of  the  nor¬ 
thern  hemisphere  of  the  core's  outer  surface.  The  presence  of  slope  or 
dip  on  the  core-mantle  boundary  would  presuppose  the  existence  of  signifi- 
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cant  variations  in  core  radius,  unless  the  interface  consisted  of  a  large 
number  of  angular  reflecting  blocks  of  relatively  small  dimensions. 


Conversion  of  Time  Residuals  to  Variations  in  Outer  Core  Radius 


In  consideration  of  the  ellipticity  correction,  equation  (6),  the 
term  due  to  Sr  at  the  core-mantle  boundary  was  deliberately  omitted.  The 
travel-time  residuals  will  now  be  equated  to  this  term  and  expressed  in 
terms  of  variation  in  core  radius.  The  proportionality  factor, 

V. 

- - 1 — : —  ,  is  plotted  as  a  function  of  epicentral  distance  in  Figure  14. 

(i 


The  angle  of  incidence  at  the  core,  ia»is  determined  by  the  relationship 


fit 


fit 


sin  ia  *  —  where  is  derived  from  the  Taggart-Engdahl  PcP  travel 
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a 

time  tables,  r  is  the  core  radius  (3477  km),  and  V  is  the  velocity 
(13.67  km/sec)  in  the  mantle  just  above  the  core.  These  values  are  from 
model  4  of  Taggart  and  Engdahl  [1968].  The  inverse  delay  factor  of 

6r 

Figure  14,  ,  increases  exponentially  with  distance  and  would  be  infinite 


at  grazing  incidence.  The  decision  to  restrict  the  fir  interpretation  to 
epicentral  distances  less  than  70  degrees  was  detc  rmined  in  part  by  this 
consideration  and  in  part  by  the  fact  that  some  of  the  better  PcP  arri¬ 
vals  seem  to  be  found  near  this  distance.  This  is  illustrated  in  the  bot¬ 
tom  half  of  Figure  13  where  PcP  from  the  Sahara  event  recorded  at  State 
College,  Pennsylvania,  at  a  distance  of  70  degrees  is  shown. 


Distribution  of  PcP  Reflection  Points 

A  computer  program  was  written  to  compute  the  point  of  reflection  for 
PcP,  assuming  that  reflection  occurs  at  the  projection  on  the  core  of  the 
midpoint  of  the  geocentric  great  circle  path  from  epicenter  to  station. 

The  program  consists  of  a  standard  geocentric  distance  and  azimuth  algo¬ 
rithm  modified  to  accept  a  subroutine  which  calculates  longitude  and  geo¬ 
centric  latitude  of  reflection  point,  azimuth  of  the  ray  at  the  reflection 
point,  the  sum  (HSUM)  of  firepicenter  Plus  6rstation  due  t0  the  ellipti¬ 
city  of  the  earth's  surface,  the  harmonic  S£  defined  in  conjunction  with 
equation  (9)  as  a  function  of  the  latitude  of  the  reflection  point,  and 
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the  sum  Sq  +  y  involves  S£  harmonics  at  the  epicenter,  re¬ 

flection  point  and  station  and  is  used  in  the  calculation  of  the  correc¬ 
tion  for  the  ellipticities  of  internal  strata  of  equal  velocity.  A  print¬ 
out  displaying  the  output  of  the  Fortran  program  as  modified  is  included 
in  Appendix  2. 

PcP  reflection  points  corresponding  to  data  used  for  determination 
of  <5rcore  are  plotted  on  an  azimuthal  equidistant  polar  projection,  and 
is  shown  in  Figure  15.  Coverage  of  the  northern  hemisphere  is  by  no  means 
complete,  but  should  be  sufficient  to  detect  broad  trends.  Coverage  of 
the  north  polar  region  is  better  than  that  near  the  equator,  and  the  ma¬ 
jority  of  the  reflection  points  fall  at  mid-latitudes. 

Presentation  of  fir  Data 

The  density  of  reflection  points  in  some  regions  precludes  the  in¬ 
clusion  of  each  point  on  a  global  residual  map.  Consequently  the  data 
are  represented  by  sector  means,  all  values  being  averaged  within  zones 
10°  of  latitude  by  10°  of  longitude.  This  method  allows  the  minimum 
dilution  of  a  large  percentage  of  the  data,  since  the  smaller  areas  en¬ 
closed  by  zones  at  higher  latitudes  are  compensated  by  the  greater  den¬ 
sity  of  reflection  points  there. 

The  sectorially  averaged  variations  in  core  radius  presented  in 
Figure  16  are  based  upon  PcP  travel  times  at  distances  less  than  70°,  cor¬ 
rected  for  elevation  and  ellipticity  only.  Averages  determined  in  the 
same  manner  for  residuals  based  upon  PcP  travel  times  with  station  and 
source  corrections  applied  are  shown  in  Figure  17.  A  similar  presenta¬ 
tion  for  fir  based  upon  PcP-P  times  is  shown  in  Figure  18. 

To  show  some  of  the  individual  data  points,  the  fir  based  upon  PcP-P 
at  reflection  points  north  of  40°  north  latitude  have  been  plotted  on  the 
larger  scale  polar  projection  of  Figure  19.  Note  the  general  agreement  of 
fir  at  neighboring  reflection  points  from  different  events. 

The  agreement  of  the  values  in  Figures  17  and  18  is  not  as  great  as 
had  been  anticipated.  As  indicated  earlier,  the  two  data  sets  should  be 
equivalent,  since  the  crust  and  uppei.  mantle  travel-time  anomalies  have, 
hopefully,  been  largely  compensated  in  each  case.  However,  since  the  error 
dependence  was  shown  to  be  different  for  the  two  sets  of  data,  the  best 
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30°  20°  10°  0  -10°  -20°  -30° 

FIGURE  15.  DISTRIBUTION  OF  PcP  REFLECTION  POINTS,  A  <  70° . 
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FIGURE  16.  SECTORIALLY  AVERAGED  6r  OF  CORE  BASED  UPON  PcP  TRAVEL 
TIMES. 
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FIGURE  19.  6r  BASED  UPON  PcP-P  DATA,  kO°  -  90°  NORTH  LATITUDE 


estimate  of  core  radius  is  obtained  by  averaging  the  sets,  producing  a  new 
set,  6r.  The  6r  was  obtained  at  each  reflection  point  by  averaging  the 
6r  determined  from  PcP  times  corrected  for  station  and  source  travel-time 
anomalies  and  the  6r  obtained  from  the  PcP-P  times,  except  in  the  few 
cases  where  station  time  was  not  known  or  no  P  arrival  time  was  available. 
Further  data  analysis  will  utilize  primarily  the  6r  set,  sectoral  averages 
of  which  are  displayed  in  Figure  20. 

Statistical  Distribution  of  the  6r  Data 

A  histogram  (Figure  21)  of  the  6r  data  was  prepared,  showing  the  num¬ 
ber  of  data  points  in  each  one  kilometer  interval.  The  data  mean  is  +0.75 
kilometers,  with  a  standard  deviation  of  20.58  kilometers.  A  few  scat¬ 
tered  data  points  more  than  two  standard  deviations  from  the  mean  fall  be¬ 
yond  the  scale  limits  of  the  abscissa  of  Figure  21  and  hence  are  not  shown. 

The  distribution  is  nearly  Gaussian,  with  little  evidence  of  bi-  or  multi¬ 
modality  which  might  indicate  a  layered  structure  consisting  of  more  than 
one  reflector.  The  distribution  of  6r  is,  ef  course,  dependent  upon  the 
spatial  distribution  of  the  reflection  points,  especially  with  regard  to 
latitude.  Any  layered  structure  of  the  mantle-core  boundary  which 
might  exist  would  be  quite  "blurred"  by  the  variation  of  core  radius  with 
latitude.  The  value  of  the  mean  of  6r  reflects  the  preponderance  of  data 
on  the  flattened  portion  of  the  core  spheroid,  north  of  35°  north  latitude. 

It  will  be  shown  that  the  data  indicate  a  somewhat  larger  increase  in  the 
core  radius  of  Taggart  and  Engdahl’s  [1968]  model  4. 

Harmonic  Representation  of  the  6r  Data 

The  distribution  and  scatter  of  the  data  of  Figures  16  through  20 
makes  contouring  difficult,  if  not  impossible.  A  method  of  smoothing  and 
interpolation  is  needed  to  present  the  data  as  a  function  of  latitude  and 
longitude,  to  extract  significant  trends,  and  to  describe  the  gross  shape 
of  the  core.  This  has  been  achieved  by  expressing  the  6r  in  a  series  of 
spherical  harmonics. 

Such  an  arbitrary  function  of  colatitude  (<j>)  and  longitude  (9)  on  the  sur¬ 
face  of  the  unit  sphere  can  be  represented  by  a  series  of  surface  harmonics, 
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FIGURE  20.  S ECTOR IALLY  AVERAGED  8r  OF  CORE. 
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11(1,0) 


■  rAjjQ  n  I 

— —  Pn(cos  40  +  l  (Anm  cos  m0  +  Bnm  sin  m6)  Pnm(cos  <j>)  (10) 

n=0  m=l 


whose  coefficients  are  given  by 


2n+l  (n-ra) J 
'Vim  _  "  2 n  (n+m) ' 


and 

_  2n+l  (n-m) ! 
nm  2 tt  (n+m) ! 


Pn  (cos  40  cos  iti0  •  sin4>  dcf>  d0 


m 

Pn  (cos  40  sin  m0  sin  <p  dip  dO 


m 

where  Pn  denotes  the  Legendre  polynomial  and  Pn  the  "associated  Lengendre 
polynomial"  defined  by 


m 

pn  00 


dm[Pn(y)] 

- jjj -  •  sinm4>,  where  y  =  cos  <f>. 

dp 


Thus  the  constants  Anm  and  Bnm  can  be  determined  by  integration  over 
the  surface  of  the  sphere.  An  alternate  approach  is  to  determine  the  A^ 
and  Bjjh,  by  a  least  squares  solution  to  a  system  of  linear  equations,  each 
of  which  has  the  form  (10).  Equation  (10)  can  be  rewritten 


H(<M> 


n 


1  1  (Vrm  Vim  + 

n=0  m=0 


nm  Vm^ 


(12) 


m 

where  Cnm  =  Pn  (y)  •  cos(m0) 
m 

and  Snm  =  Pn  (y)  •  sin(m0) 


(13) 


are  the  so  called  "tesseral"  surface  harmonics  and  Anm  and  Bnm  are  suit¬ 
ably  chosen  constants.  The  Vim  and  Snm  have  been  normalized  such  that 


r  2  TT  -7T  .2TT  -TT 

J  j  c^m  sin  <t>  d(|)  d0  =  J  j  S^m  sin  4>  d0  d0  =»  Att  . 

The  normalizing  factor,  Fnm  is vj^n+^'  for  the  zonal  harmonics  (m*0)  and 
v/ 2 (2n+l) (nam) I /(n+m) !  for  the  tesseral  and  sectoral  harmonics. 

Each  data  point  defines  values  of  <J>,  0,  and  H  and  thus  determines  an 
equation  of  the  form  (12).  The  6r  data  were  expressed  as  a  series  of 
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surface  harmonics  through  degree  (N)  five.  The  algorithm  used  is  a  modi¬ 
fication  and  extension  of  Pollack's  [1968]  program  utilizing  MacMillans 
[1930]  tabulation  of  the  tesseral  harmonics  through  degree  four  in  Car¬ 
tesian  coordinates.  The  program  as  modified  can  be  used  to  determine 
tesseral  harmonic  coefficients  through  degree  and  order  five  and  zonal 
harmonic  coefficients  to  higher  degrees.  It  also  incorporates  another 
technique,  described  herein  as  "step  harmonics",  found  to  be  especially 
useful  in  the  analysis  of  data  which  are  not  uniformly  distributed  on  the 
sphere.  The  lack  of  sufficient  constraints  in  the  form  of  data  in  the 
equatorial  regions  and  southern  hemisphere  proved  to  be  a  serious  problem 
in  the  harmonic  analysis.  The  least  squares  "best  fit"  to  the  observa¬ 
tions  obtained  in  the  initial  attempts  resulted  in  some  of  the  coefficients 
becoming  very  large.  While  the  4th  or  5th  degree,  series  thus  generated 
provided  a  slightly  better  fit  to  the  data  than  was  later  obtained  by  the 
step  harmonic  technique,  the  results  were  unacceptable  because  they  pro¬ 
duced  a  core  which  protruded  from  the  earth's  surface  near  the  south  pole. 
Constraints  were  imposed  in  the  form  of  artificial  data  points,  6r  =  0, 
introduced  at  strategic  locations  in  the  southern  hemisphere,  but  this 
was  found  to  be  both  mathematically  and  psychologically  unsuitable.  The 
method  of  step  harmonics  consists  of  the  independent  determination  of  co¬ 
efficients,  a  few  terms  at  a  time.  The  contribution  of  these  terms  is  re¬ 
moved  from  the  data  and  the  residual  used  to  compute  the  next  set  of  co¬ 
efficients.  The  coefficients  computed  by  this  method  will  depend  upon 
the  order  of  evaluation  and  thus  will  not  be  unique  although  the  resulting 
representations  of  the  <5r,  within  the  limits  of  the  data  scatter  and  dis¬ 
tribution  of  reflection  points,  will  be  similar.  Since  any  least  squares 
technique  is  highly  sensitive  to  extreme  data  values,  the  6r  data  set  was 
truncated  at  ±20  km  (approximately  one  standard  deviation) ,  reducing  the 
data  points  to  223.  First  degree  terms  were  omitted  in  the  analysis. 

This  is  equivalent  to  the  assumption  that  the  surface  of  the  core  and  the 
earth's  surface  are  concentric.  Two  sets  of  harmonic  coefficients  were 
computed  using  different  step  sequences.  These  values  are  shown  in  Table 
4.  Set  1  was  determined  by  first  evaluating  Cqq  and  the  second  degree 
terms.  Subsequent  steps  determined  the  third,  fourth,  and  finally  the 
fifth  degree  coefficients.  The  first  step  in  the  determination  of  co¬ 
efficient  set  2  was  the  computation  of  the  zonal  harmonic  coefficients 
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Table  4.  Harmonic  Coefficients 


Coefficient 


Set  1 


Set  2 


A 

A 

A 

B 

A 

B 

A 

A 

B 

A 

B 

A 

B 

A 

A 

B 

A 

B 

A 

B 

A 

B 

A 

A 

B 

A 

B 

A 

B 

A 

B 

A 

B 


00 

20 

21 

21 

22 

22 

30 

31 

31 

32 

32 

33 
33 

40 

41 

41 

42 

42 

43 

43 

44 
44 

50 

51 

51 

52 

52 

53 

53 

54 

54 

55 
55 


1.902 

-1.822 

-1.220 

1.180 

-2.186 

-0.676 

-0.137 

0.499 

0.020 

0.257 

0.770 

-0.740 

0.155 

-0.521 

0.559 

-0.064 

-0.131 

0.366 

-1.163 

-0.286 

0.275 

0.400 

-0.481 

0.210 

-0.295 

-0.502 

0.098 

-0.249 

-0.063 

0.809 

0.554 

0.656 

-0.164 


1.229 

-1.499 

-0.237 

0.122 

-0.213 

-0.349 

1.338 

-0.214 

0.266 

-0.005 

-0.076 

0.429 

0.400 

-1.410 

-0.185 

0.716 

0.093 

0.182 

-0.045 

0.407 

-1.029 

-0.492 

-0.365 

0.797 

-0.376 

-1.154 

-0.019 

-1.347 

-0.429 

1.359 

0.550 

1.614 

-1.348 
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^00’  ^20’  ^30*  anC^  \o*  Their  contribution  was  then  subtracted  from  the 
6r  data  and  the  resulting  residuals  used  in  the  determination  of  the  5th 
degree  terms.  This  procedure  was  repeated  to  compute  the  remaining  fourth, 
third,  and  second  degree  terms.  The  series  of  spherical  harmonics  cor¬ 
responding  to  sets  1  and  2  describe  the  existing  data  without  introducing 
profound  anomalies  in  regions  where  there  are  no  data.  It  must  be  kept  in 
mind  that  the  principal  reason  for  the  harmonic  expansion  was  the  smooth¬ 
ing  and  interpolation  of  the  data.  The  signs  of  and  the  second  de¬ 
gree  coefficients  in  Table  4  are  in  agreement  and  may  be  representative  of 
the  gross  size  and  shape  of  the  core.  The  other  terms  are  included  only 
to  improve  the  representation  of  the  data.  The  rms  error  in  this  repre¬ 
sentation  is  approximately  7.5  kilometers  for  each  set.  The  convention 
of  positive  west  longitude  used  in  the  figures  throughout  this  paper  was 
reversed  in  the  computation  of  the  harmonic  coefficients,  rendering  east 
longitude  positive  in  keeping  with  standard  practice. 

An  Estimate  of  the  Shape  of  the  Core 

The  values  of  the  terms  A£q  and  computed  in  the  first  steps  of 

the  analysis  determining  coefficient  sets  1  and  2  show  better  agreement 
between  sets  than  the  other  terms.  Although  obvious  from  rotational  con¬ 
siderations,  the  oblateness  of  the  mantle-core  boundary  has  hitherto  not 
been  demonstrated  on  the  basis  of  PcP  travel  times.  The  ellipticity  of 
the  core  was  determined  from  the  A2Q  coefficient  by  evaluating  the  product 
^20  ^20  at  e9uator  anc*  P°le*  The  ellipticities  obtained  are  .00176  for 
coefficient  set  1  and  .00145  for  set  2.  These  values  are  notably  smaller 
than  the  .00260  calculated  by  Bullen  [1936]  on  the  basis  of  P-  and  S-wave 
velocities,  using  the  Williams on- Adams  method  and  the  Radau- Darwin  approxi¬ 
mation.  Additional  PcP  observations  are  required  in  the  equatorial  re¬ 
gions  to  confirm  this  discrepancy. 

The  Aqq  coefficients  in  Table  4  correspond  to  an  increase  of  the 
mean  core  radius  of  the  Taggart-Engdahl  Model  4  of  about  1.5  kilometers. 

The  6r  as  defined  by  the  coefficient  sets  1  and  2  have  been  contour¬ 
ed  and  are  shown  in  Figures  22  and  23  respectively.  The  contours  overlie 
Figure  15  to  indicate  the  areas  where  data  exist.  The  variation  in  the 
6r  is  about  20  kilometers  in  each  case,  and  the  anomaly  patterns  are  very 
similar.  It  will  be  shown,  however,  that  the  potential  fields  at  the 
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surf. ico  resulting  from  coefficient  sets  1  and  2  are  quite  different.  The 
r  distributions  of  Figures  22  and  23  bear  little  resemblance  to  those 
determined  by  Vogel  [1%0|  using  core  reflections  from  earthquake  sources. 

Computation  of  Geoid  Heights 

Because  of  the  large  density  jump  believed  to  exist  at  the  mantle- 
core  boundary,  broad  undulations  on  the  order  of  kilometers  should  result 
in  observable  gravity  variations  at  the  earth's  surface.  The  correspond¬ 
ing  geo  id  heights  may  be  determined  by  considering  a  surface  density, q,  at 
the  mantle-core  boundary,  defined  by 

q(<MO  =  Ap  •  6r  (<J>,G) 


where  Ap  is  the  density  contrast  at  the  mantle-core  boundary.  This  may 
be  written  in  terms  of  the  harmonic  expansion  of  <$r, 

oo  n 

q  =  Ap  y  y  (A_m  C  +  B  S__)  . 
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Following  Jeffreys  [1962],  the  gravitational  potential  anomaly  at  the 
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where  r  and  R  are,  respectively,  the  radii  of  the  core  and  of  the  earth,  and 
G  is  the  universal  gravitational  constant.  The  geoid  height  is  then 

Ah  (4> ,  0)  =  AW^-9-^  . 

The  hypothetical  geoid  heights  resulting  from  the  6r  of  Figure  22  repre¬ 
sented  by  coefficient  set  1,  Table  4,  omitting  the  C^q  term  and  assuming 
a  value  of  A  gm/cc  for  Ap,  are  contoured  in  Figure  24.  The  sectoral  har¬ 
monic  S dominates  a  height  distribution  characterized  by  variations  of 
*300  meters.  The  geoid  heights  corresponding  to  the  6r  of  Figure  23,  com¬ 
puted  from  coefficient  set  2,  are  shown  in  Figure  25.  The  anomalies  are 
more  local  in  extent  and  smaller  in  amplitude  than  those  of  Figure  24. 
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FIGURE  2b.  GEOID  HEIGHTS  COMPUTED  FROM  Br  DISTRIBUTION  OF  FIGURE  22. 
(Contour  interval  =  40  meters) 
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Geoid  heights  computed  by  Kaula  [1966]  using  harmonic  coefficients 
through  the  8th  degree  derived  from  observations  of  satellites  are  shown 
in  Figure  26.  The  contours  were  transferred  to  the  polar  projection  to 
facilitate  comparison  with  Figures  24  and  25.  While  the  range  of  values 
in  Figures  25  and  26  is  about  the  same,  their  spatial  correlation  is  poor. 
Conspicuous  by  its  absence  in  Figure  24  is  the  positive  geoid  anomaly  cen¬ 
tered  on  the  Greenwich  meridian  at  50°N  latitude.  The  lack  of  agreement 
could  be  due  to  many  factors.  Among  these  are: 

1.  1.  The  travel-time  variations  used  in  computing  the  5r  may  be 

the  result  of  velocity  variations  in  the  mantle  rather  than  changes  in 
depth  to  the  mantle-core  boundary. 

2.  The  inferred  irregularities  of  mass  distribution  at  the 

mantle-core  boundary  may  be  masked  by  extensive  irregularities  at  other 
depths.  Total  or  partial  compensation  may  occur.  Such  compensation  would 
have  to  o„cur  with  regard  to  the  contribution  of  the  second  degree  sec¬ 
torial  harmonics an<^  ^21’  Jeffreys  [1962]  has  shown  must  be  ab¬ 

sent  in  the  potential  since  Z  is  the  axis  of  greatest  moment  of  inertia  of 
the  earth. 

3.  The  density  contrast  at  the  mantle-core  boundary  may  be 
less  than  the  value  of  4  gm/cc  used  in  the  geoid  height  computations. 
Buchbinder  [1968b]  has  built  a  case  for  continuity  of  density  across  the 
boundary  on  the  basis  of  an  amplitude  minimum  and  phase  reversal  of  PcP 
at  epicentral  distances  near  32°. 

4.  The  velocity  at  the  base  of  the  mantle  may  be  lower  than 
that  estimated  by  Taggart  and  Engdahl  [1968],  Decreasing  this  veloci¬ 
ty  would  tend  to  decrease  6r  and  thus  decrease  the  hypothetical  geoid 
heights.  Of  the  84  surface  focus  PcP  arrival  times  used  by  Taggart  and 
Engdahl,  only  four  were  at  epicentral  distances  greater  than  65°.  Al¬ 
though  PcP  travel  times  at  shorter  distances  are  not  very  sensitive  to 
lower  mantle  structure,  analysis  of  these  data  indicated  a  slight  pre¬ 
ference  for  a  model  exhibiting  a  gradual  increase  of  velocity  with 
depth  in  this  region.  Taggart  and  Engdahl' s  [1968]  PcP  travel  time 
tables  are  based  on  this  model  with  a  core  radius  of  3477  km.  The  tra¬ 
vel  times  tabulated  in  Appendix  1,  some  of  which  are  taken  from  their 


45 


30°  20°  10°  0  -10°  -20°  -30° 


FI  JURE  ?6.  GEO ID  HEIGHTS  COMPUTED  FROM  SATELLITE  OBSERVATIONS  (AFTER 

KAULA  [1966]). 

(Contour  interval  =  20  meters) 


46 


A  ■  96°  in  terms  of  a  core  radius  of  3550  km.  This  interpretation  is  not 
consistent  with  PcP  travel  times,  which  are  sensitive  to  major  changes  in 
core  radius.  It  would  appear  that  this  diffraction  boundary  is  associated 
with  a  decrease  of  velocity  with  depth  in  the  lowest  100  km  of  the  mantle, 
rather  than  with  the  position  of  the  mantle-core  boundary.  Ergin  [1967] 
has  found  such  a  velocity  decrease  in  the  lower  mantle  necessary  to  ex¬ 
plain  the  range  of  observations  for  PKP2  phases. 

The  sharpness  of  the  shadow  boundary  at  periods  near  one  second,  as 
determined  by  Sacks  [1967]  using  earthquake  data  and  as  indicated  by  the 
data  of  Figure  28,  implies  that  the  amplitudes  of  these  short-period 
(<1  sec)  diffracted  P  phases  are  probably  too  small  beyond  about  100°  to 
be  detected. 

Such  short-period  arrivals  do  exist,  however,  beyond  100°  and  have 
been  reported  to  distances  of  134®  or  more  (Ergin,  1967,  Lehmann,  1953). 

Bolt  [1969]  has  observed  that  first  arrivals,  designated  "diffracted  PcP", 
in  the  distance  range  of  100°  to  115°  fall  on  a  line  of  slope  dt/dA  of 
4.60  seconds  per  degree.  His  lower  mantle  P-wave  velocities,  based  on 
this  slope  and  on  PcP  travel  times,  decrease  to  a  value  of  13.3  km/ sec 
at  the  mantle-core  boundary.  Additional  evidence  of  a  velocity  decrease 
at  the  base  of  the  mantle  may  be  determined  from  the  "cutoff  puxnts"  in¬ 
dicated  by  the  short-period  core  phase  travel-time  data  of  Engdahl  [1968]. 

The  cutoff  distance  of  PcP  must  be  equal  to  twice  the  cutoff  distance  of 
PKKP  reduced  by  the  cutoff  distance  of  PKKKKP.  This  implies  PcP  propa¬ 
gation  to  a  distance  of  124°. 

Presentation  of  Amplitude  Data 

The  amplitude/period  (A/T)  ratio  of  P  and  PcP  phases  and  the  A/T  ratio 
PcP/P  are  tabulated  in  Appendix  1  for  a  number  of  arrivals  used  in  the  analy¬ 
sis  of  PcP  and  P  travel  times.  The  variation  of  PcP  A/T  (vertical  component) 
ratios  with  distance  is  shown  in  Figure  27.  All  events  have  been  normalired  to 
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can  only  be  estimated  since  the  lower  amplitude  arrivals  are  either  lost 
in  the  P  coda  or  are  below  the  threshold  of  detectability.  The  range  of 
scatter  is  probably  on  the  same  order  as  that  of  P  for  comparable  paths. 

The  vertical  amplitude/period  ratio  for  P  waves  from  nuclear  explosions 
in  the  distance  range  70°  to  105°  ip  shown  in  Figure  28,  normalized  to 
magnitude  6.5.  Although  the  variation  exceeds  40  decibels,  the  decrease 
of  amplitude  at  a  distance  of  about  96°  is  apparent,  marking  the  beginning 
of  the  "shadow  zone"  for  short-period  P  waves.  The  ratio  PcP/P  (vertical 
components)  shown  in  Figure  29  increases  somewhat  with  distance  as  indi¬ 
cated  by  the  dashed  logarithmic  mean.  The  period  ratio  of  P  to  PcP  is 
shown  in  Figure  30  as  a  function  of  epicentral  distance.  The  predominance 
of  values  less  than  unity  for  the  Novaya  Zemlya  ratios  is  probably  the 
nost  significant  feature  of  this  plot.  The  PcP  amplitude/period  ratio  is 
presented  as  a  function  of  6r  in  Figure  31.  The  relationship  appears  to 
be  a  random  one,  as  does  that  of  the  PcP/P  ratio  A/T  as  a  function  of  the 
PcP-P  travel-time  residual  shown  in  Figure  32. 

A  Model  Study  of  Crust  and  Upper  Mantle  Effects 

The  high  degree  of  scatter  in  the  PcP  amplitude  and  travel-time  data 
tends  to  obscure  relationships  due  to  conditions  at  or  near  the  mantle-core 
boundary.  Although  lateral  variations  in  velocity  may  exist  in  the  lower 
mantle  (Chinnery  [1969])  it  is  probable  that  most  of  the  observed  data 
scatter  is  due  to  the  inhomogeneity  of  the  crust  and  upper  mantle.  Broad 
regional  variations  in  structure  produce  travel-time  anomalies  which  may 
be  compensated  by  the  use  of  source  and  station  corrections  or  the  PcP-P 
time  difference.  Anomalies  due.  to  more  abrupt  changes  in  structure  are 
difficult  to  correct,  requiring  detailed  knowledge  of  crust  and  upper  man¬ 
tle  structure  or  empirical  determination  of  corrections  as  functions  of 
azimuth  and  epicentral  distance.  The  velocity  structure  shown  in  Figures 
33  and  34  is  similar  to  the  upper  mantle  low  velocity  zone  proposed  by 
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FIGURE  51.  PcP  A/T  RATIO  AS  A  FUNCTION  OF  6r. 
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FIGURE  32.  PcP/P  RESIDUAL  OF  A/T  AS  A  FUNCTION  OF  THE  PcP-P  TRAVEL 
TIME  RESIDUAL. 
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FIGURE  33.  EFFECT  OF  VARIATIONS  IN  CRUST  AND  UPPER  MANTLE  STRUCTURE 
ON  THE  PcP-P  TRAVEL  TIME  DIFFERENCE. 
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DEPTH  IN  KILOMETERS  VARIATION  OF  AMPLITUDE  RATIO  PcP/P 


Nuttli  and  Bolt  [1969]  to  account  for  P-wave  travel  residuals  at 
northern  California  stations.  The  lateral  variations  are  probably  the 
most  extreme  proposed  anywhere  with  the  exception  of  a  few  island  arc 
areas.  To  evaluate  the  effect  of  such  a  model  on  the  PcP-P  time  interval, 

PcP  and  P  were  propagated  through  the  structure  using  the  digital  ray 
tracing  program  developed  by  Jackson  [1968].  Plane  waves  and  a  flat  earth 
were  assumed,  with  a  10°  angle  of  incidence  at  the  base  of  the  low 
velocity  zone  assumed  for  PcP  (Figure  33),  and  35°  for  P  (Figure  34).  These 
are  approximately  the  appropriate  angles  at  an  epicentral  distance  of  35 
degrees.  Travel-time  curves  for  PcP  and  P  were  computed  from  the  ray 
travel  times.  The  changes  in  the  PcP-P  time  interval  due  to  lateral  vari¬ 
ations  of  the  velocity  structure  are  shown  at  the  top  of  Figure  33.  It 
is  evident  from  an  examination  of  Figures  33  and  34  that  there  are  zones 
of  energy  concentration  and  shadow  zones.  The  travel-time  differences 
of  the  overlapping  wave  fronts  are  generally  a  small  fraction  of  the 
period  at  1  Hz,  so  the  arrivals  will  be  nearly  in  phase.  The  PcP  and  P 
vertical  amplitudes  were  taken  to  be  proportional  to  the  ray  density  at 
the  model  is  shown  at  the  top  of  Figure  34.  A  window  of  48  km  was  used 
to  determine  the  ray  density,  thus  smoothing  the  amplitude  variations  approx¬ 
imate  diffraction.  These  variations,  which  were  found  to  be  as  large  as 
12  db  for  P  and  8  db  for  PcP,  could  thus  result  in  a  variation  of  as  much 
as  20  db  for  the  ratio  PcP/P.  A  similar  crust  and  upper  mantle  structure 
to  the  north  in  British  Columbia  could  account  in  part  for  the  anomalous 
amplitudes  of  P  and  PcP  observed  for  the  Longshot  event  at  Prince  George, 
British  Columbia,  Figure  6.  However,  it  is  likely  that  the  structure  near 
Longshot  could  also  contribute  to  the  amplitude  anomaly. 

Summary  and  Conclusions 

Travel  times  of  PcP  and  P  phases  from  nuclear  and  high  explosive 
sources  have  been  interpreted  in  terms  of  variations  in  the  radius  of  the 
outer  core.  This  interpretation  results  in  a  core  which  is  slightly  lar¬ 
ger  than  that  of  the  reference  model  of  Taggart  and  Engdahl  [1 968 J  and  has 
less  ellipticity  than  estimated  by  Bullen  [1936],  The  distribution  of  data 
does  not  permit  a  conventional  spherical  harmonic  analysis  of  the 
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inferred  variation  in  core  radius.  A  modified  spherical  harmonic  analy¬ 
sis  method  is  devised  to  smooth  the  data  and  to  estimate  the  shape  of  the 
core.  The  variation  of  core  radius  determined  from  this  representation  is 
approximately  f10  kilometers  on  the  basis  of  terms  to  degree  and  order  5. 

The  standard  deviation  of  the  data  is  approximately  ±20  kilometers,  in¬ 
dicating  the  presence  of  variations  of  higher  degree  and  order,  many  of 
which  are  probably  due  to  factors  other  than  undulation  of  the  mantle-core 
boundary.  Geoid  heights  computed  from  the  harmonic  coefficients  (Set  2)  do 
not  correlate  with,  but  are  within  the  range  of,  those  determined  from  gravity 
data.  Analysis  of  PcP  amplitude  data  shows  no  significant  variation  of 
amplitude  with  distance  and  an  amplitude  scatter  on  the  same  order  as  that 
of  P.  The  possible  contribution  of  an  upper  mantle  low  velocity  zone  of 
varying  thickness  to  the  amplitude  and  travel-time  scatter  is  demonstrated. 

Recommendations 

The  present  study  was  greatly  limited  by  the  geographical  distribution 
of  events.  Contained  nuclear  explosions  and  underwater  high  explosive 
shots,  especially  those  whose  zero  times  and  locations  are  made  known,  pro¬ 
vide  a  unique  seismo logical  tool  for  the  determination  of  the  detailed 
structure  of  the  earth.  One  of  the  peaceful  uses  of  atomic  energy  should 
be  the  detonation  of  underground  explosions  of  approximate  magnitude  6  at 
locations  selected  to  provide  eventual  global  coverage.  These  detonations 
should  be  announced  well  in  advance  to  allow  maximum  participation  in  the 
recording  program.  Shots  in  the  central  Pacific  and  in  the  southern 
hemisphere  are  especially  needed  to  fill  the  largest  existing  data  gaps. 

Additional  surface- focus  sources  will  permit  refinement  of  source  and 
station  corrections  as  functions  of  both  azimuth  and  epicentral  distance. 

The  period  ratio  of  P  to  PcP  was  found  to  be  generally  less  than  one 
for  the  Novaya  Zemlya  event  of  October  27,  1966.  For  other  surface- focus 
events  this  ratio  has  been  predominantly  equal  to  or  greater  than  one. 

Events  in  this  region  warrant  further  investigation  to  determine  whether 
this  unusual  phenomenon  is  due  to  the  velocity  structure  beneath  the 
source  or  to  a  complex  source  function. 

Travel  times  and  amplitudes  of  short-period  P  /aves  in  the  distance 
range  90°  to  140°  and  of  PcP  beyond  70°  should  be  studied  to  determine  in 
detail  the  velocity  structure  of  the  lowest  few  hundred  kilometers  of  the 
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mantle  and  the  magnitude  of  the  density  jump  at  the  mantle-core  boundary. 
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Event:  Semipalatinsk.  1/15/65  06:00:00.0  49. 92N  78.92E 
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Event:  Chase  IV 
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PcP  TRAVEL  TIME  REFERENCES 
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Note :  Amplitude  and  period  data  included  is  generally  from  the  same 
author  as  the  travel-time  data.  Most  of  the  Longshot  values  of  PcP 
A/T  and  the  PcP/P  ratio  A/T  are  read  from  amplitude-distance  plots  pre¬ 
pared  by  D.  S.  Carder  (personal  communication). 
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